Astronomy & Astrophysics manuscript no. 


February 2, 2008 


(DOI: will be inserted by hand later) 





Numerical simulation of the surface flow of a companion star in 

a close binary system 

Kazutaka Oka 1 , Takizo Nagae 1 , Takuya Matsuda 1 , Hidekazu Fujiwara 2 and Henri M.J. Boffin 3 



(N 
O 
O 
(N 

< 



> 
o 
o 

(N 
oo 

o 

O 

S3 

Oh: 

6 



X 



Department of Earth and Planetary Sciences, Kobe University, Kobe 657-8501, Japan 

e-mail: kazutaka@kobe-u.ac.jp, nagae@kobe-u.ac.jp, tmatsuda@kobe-u.ac.jp 

IBM Japan Ltd., Yamato-shi, Kanagawa 242-8502, Japan 

Royal Observatory of Belgium, 3 av. Circulaire, 1180 Brussels, Belgium 

e-mail: Henri.Boffin@oma.be 



Received 29 April 2002/ Accepted 12 July 2002 

Abstract. We simulate numerically the surface flow of a gas-supplying companion star in a semi-detached binary 
system. Calculations are carried out for a region including only the mass-losing star, thus not the mass accreting 
star. The equation of state is that of an ideal gas characterized by a specific heat ratio 7, and the case with 
7 = 5/3 is mainly studied. 

A system of eddies appears on the surface of the companion star: an eddy in the low pressure region near the 
LI point, one around the high pressure at the north pole, and one or two eddies around the low pressure at the 
opposite side of the LI point. Gas elements starting near the pole region rotate clockwise around the north pole 
(here the binary system rotates counter-clockwise as seen from the north pole). Because of viscosity, the gas drifts 
to the equatorial region, switches to the counter-clockwise eddy near the LI point and flows through the LI point 
to finally form the LI stream. 

The flow field in the LI region and the structure of the LI stream are also considered. 
Key words, binaries : close - star : evolution - stars : mass loss - accretion, accretion disks 



1. Introduction 

A semi-detached binary is a system consisting of a pri- 
mary star and a Roche-lobe filling secondary star. In cat- 
aclysmic variables (CVs), for instance, the primary star is 
a white dwarf, while the Roche-lobe filling star is a low- 
mass main-sequence star, from which gas is supplied to 
the primary star via an accretion disc (see Warner 1995 
for a review). 

Accretion discs have been the main target of astro- 
physical research, while gas supplying companion stars 
has not been investigated fully from both theoretical and 
observational viewpoints. Knowledge of the properties of 
the companion star is required in order to understand the 
evolutionary mechanism of the binary systems. One rea- 
son that the companions have not attracted much atten- 
tion from astrophysicists is that they had little observ- 
able c onstraints. Recently, however, Dhillon and Watson 
( 200C ) have demonstrated the technique of imaging the 
companion stars in CVs using Roche tomography. 

The companion stars in CVs are similar to low-mass 
main-sequence stars in their gross properties. However, 
companion stars exist in a different environment from 
those of isolated stars. They might be influenced by ef- 



fects like irradiation from the accretion disc around the 
primary star, its rapid rotation, Roche-lobe shape, mass 
loss, and so on. These effects could in turn affect the evo- 
lution of the binary systems. 



1.1. Astrostrophic wind of Lubow and Shu 



Lubow and Shu (1975, 1976) conducted pioneering work 
in the study of gas dynamics in semi-detached binaries. In 
their paper, they analyzed the surface layer of the compan- 



ion star, as well as the LI region (Lubow and Shu 1975) 
and the LI stream (Lubow and Shu 1976). In their analy- 
sis they solved the hydrodynamic equations by a matched 
expansion technique and reduced partial differential equa- 
tions to ordinary differential equations. Their method was 
superior to numerical methods at that time considering 
the poor capability of computers. Since then, the progress 
in computer power has been considerable, so that it is 
worthwhile to investigate the same problem numerically. 
Our main target in the present paper is to investigate the 
surface layer, the LI region and the LI stream numerically. 
The region around the accreting compact star has been in- 
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vestigated separately by our group (Makita, Miyawaki and 
Matsuda pOOC| , Matsuda et al. |2000| , Fujiwara et al. |200l| ). 

Lubow and Shu ( 1975| ) stated that "the horizontal 
component of the fluid velocity is parallel to isobars in the 
astrostrophic approximation. Like the large-scale circula- 
tion patterns in the atmosphere and oceans of the Earth, 
the flow does not simply proceed from high pressures to 
low, but is directed around the contours of equal pressure. 
In the northern hemi-lobe, the flow proceeds by keeping 
the region of high (low) pressure to the right (left); in 
the southern hemi-lobe, by keeping high (low) pressure to 
the left (right)." On a given equi-potential surface they 
expected "the pressure to be largest at the poles (for the 
same reason that meridional circulation occurs in rotating 
stars) and smallest at LI (because of the mass-loss flow)" ; 
therefore "the excess circulation should be counter to the 
sense of the orbit rotation." 

They predicted that "on each equipotential surface, 
the flow may be astrostrophic, i.e., parallel to isobars, but 
a net transfer of matter within the entire surface layer 
to the equator is possible ... if there is an outward flux 
of matter at the bottom of the surface layer." They also 
considered that "another effect which reinforces the con- 
cept that there must be a slow drift of matter from the 
poles to the equator, and thence to LI, occurs if the enve- 
lope of the contact component is convective. In this case, 
the effective fractional drag exerted on the upper layers by 
the lower layers must produce the "Ekman effect" (Ekman 



1922) where the induced secondary flow corresponds to a 



drift across isobars from high pressures to low." 

One of the aims of the present paper is to investigate 
the flow pattern of the gas in the surface layer of the com- 
panion numerically and to compare the resul ts with the 
qualitative prediction of Lubow and Shu ( 1975 ). This kind 
of problems may be termed stellar meteorology. The struc- 
ture of the LI region and the LI stream are also consid- 
ered. 



1.2. Our previous work 

Two-dimensional finite- volume hydrodynamic simulations 
of flows in a semi-detach ed bi n ary w ere done by Sawada, 
Matsuda and Hachisu ( |l986| , |l987 ). Three-dime nsiona l 
calculation was done by Sawada and Matsuda ( 1992| ). 
They were concerned mainly by the accretion disc and 
did not pay much attention to the surface flow on the 
secondary companion. 



Recently, Fujiwara et al. (2001) performed three- 
dimensional hydrodynamic simulations of a semi-detached 
binary system including the companion star. The equation 
of state is that of an ideal gas characterized by the specific 
heat ratio 7 and the case with 7 = 1.01, meaning that the 
gas is almost isothermal, was mainly studied. The reason 
for the choice of 7 = 1.01 is to mimic a radiative cooling 
effect occurring in an accretion disc. Although their main 
concern was the hydrodynamics of the accretion disc, they 
also obtained the flow pattern on the surface of the com- 



panion. Their results showed that gas flows from the pole 
region to the equatorial region rather directly, although 
small eddy patterns appeared on the surface of the com- 
panion star. One, however, has to be careful to select the 
values of the specific heat ratio. In the surface region of 
the companion, where the gas is expanding outward, the 
temperature of the gas is reduced due to the adiabatic 
expansion. Thus the case with 7 = 1.01, which results in 
the input of energy into the expanding gas, may not be 
appropriate for the present purpose. 

This paper is organized as follows. In § 2 we describe 
the numerical method and physical assumptions. In § 3 
we show the results of our numerical simulations. In § 4 
we discuss the dynamics of the obtained eddy system. A 
summary is given in § 5. 

2. Assumptions and Numerical Method 

2.1. Assumptions 

We consider the present study as a first step in the investi- 
gation of the surface flow of the companion, and therefore 
do not take into account such complex effects as viscosity 
(except numerical), magnetic fields, and irradiation from 
the accretion disc. The equation of state considered here 
is that of an ideal gas characterized by a specific heat ratio 
7, and the cases with 7 = 5/3 and 7/5 are mainly studied, 
although the cases of 7 = 1.2 and 1.01 are also considered 
to compare the present results with those of Fujiwara et 
al. ( EOOlD . 

In the present study, the binary system consists of a 
primary star with mass Mi and a companion star with 
mass M% rotating counter-clockwise when viewed from the 
north pole. The mass ratio of the two stars, q = Mi/M\, 
is assumed to be unity in the present work. We normalize 
all physical quantities as follows: the length is scaled by 
the separation A between the centres of the two stars. 
The time is scaled by l/li, where f2 denotes the orbital 
angular velocity, and therefore the orbital period becomes 
27T. The density at the inner boundary is taken to be unity. 
The gravitational constant G is eliminated using the above 
normalization. 



2.2. Numerical method 

We use a Cartesian coordinate system. The origin of the 
coordinate is located at the centre of the companion star. 
The x— axis coincides with the line joining the centre of 
the two stars. We set the x — y plane as the orbital plane, 
thus the z— axis is perpendicular to the orbital plane and 
is oriented in the same direction as the angular momentum 
vector of the orbital rotation. The LI point is located at 
(0.5, 0, 0) and the primary star is located at (1, 0, 0), which 
is out of the present computational region. 

The computational region is a rectangular box of size 
-0.5 < x < 0.7, -0.5 < y < 0.5 and < z < 0.5. Note 
however that part of this box contains the mass-losing 
star, for which we assume that the inner part contains 
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uniform density gas at rest (see Sect. 2.3). With the as- 
sumption of symmetry of physical quantities around the 
orbital plane, the calculations are performed only in the 
upper half region above this plane. This region is divided 
in 121x101x51 grid points. In order to test the assump- 
tion of symmetry about the orbital plane, we performed 
a test in which the whole space was taken into account 
in the computation. We confirmed that the assumption of 
the symmetry was good enough. 

In the present study, we use t he sa me scheme as in 
Makita, Miyawaki, and Matsuda (200C), Matsuda et al. 
( 200C ) and Fujiwara et al. ( 2001 ): the simplified flux split- 
ting (SFS) finite- volume method proposed by Jyounouchi 
et al. ( 1993 ) and Shima and Jyouno uchi (1994 ). We refer 
to Makita, Miyawaki, and Matsuda (200C) for the numer- 
ical and the SFS schemes. With this MUSCL-type tech- 
nique, we can keep the spatial and temporal accuracy at 
second-order levels. 



2.3. Boundary and initial conditions 

We apply the same boundary conditions as in Fujiwara 
et al. (2001) except for the inner boundary location of 
the companion star. The inner boundary is assumed to be 
an equi-potential surface slightly smaller than the critical 
Roche lobe: the mean radius of the inner boundary is cho- 
sen to be 0.45 in our dimensionless unit. Note that this 
inner boundary does not necessarily mean the real sur- 
face of the companion, but is chosen only for numerical 
purposes. We have also made simulations where the inner 
boundary was chosen to be 0.40 and 0.50 and the results 
are very similar. In the rest of this paper, we will thus 
concentrate on the case with an inner boundary of 0.45. 



What we call here the "surface!' of the companion star 
might not correspond to the photosphere, where optical 
depth reaches unity. Lubow and Shu ( p. 975 ) pointed out 
that "the depth over which the optical depth at visual 
wavelengths reaches unity is not a priori related to the 
depth over which appreciable mass flow occurs." In the 
present study, we are interested in the surface flow of the 
companion, and therefore we define the "surface" as the 
region where appreciable flow occurs. 

The inside of the companion star is filled with gas hav- 
ing zero velocity, density po = 1, and sound speed c Q , 
where c is a free parameter. We assume cq = 0.05 in the 
present work, although a case with Co = 0.1 is also cal- 
culated and does not show much difference in the results. 
Such a value for the sound speed might, at first glance, 
look very large, if one would only take into account the ef- 
fective temperature of the star and scale the corresponding 
speed with the orbital velocity. This is however not correct 
because as we noted above, the inner boundary is not the 
"surface" of the star but may well lie more inside the star, 
where the gas is much hotter. Temperatures of several 10 5 
K corresponding to sound speeds of the order of 50 km/s 
or more are therefore not unrealistic. 



The assumption that the gas inside the companion star 
has zero velocity does not mean that there is no gas out- 
flow from the companion star. The gas does flow out when 
the pressure above the inner boundary is lower than that 
inside the companion. The outflow flux is calculated by 
solving Riemann problems at each step. 

The outside of the outer boundary is assumed to be 
always filled with a gas with velocity 0, density pi, and 
sound speed ci, where p\ and c\ are parameters. Use of 
these boundary conditions ensures a stable calculation. 
Note that inflow or outflow of the gas through the outer 
boundary is also possible, and the flux is calculated by 
solving Riemann problems between the inside and outside 
of the computational region at each step. 

At the initial time t = 0, the entire region, except the 
inside of the companion, is occupied by a gas with velocity 
0, density pi — 10~ 5 , and sound speed c\ = ^/l0. 

To know how gas elements from the companion build 
the LI stream is interesting, because this stream does 
bring the material from the companion to the primary 
via the LI point and an accretion disc. Thus we investi- 
gate the structure of the LI stream as well as the surface 
flow in the present work. 

Simulations are run up to t = 62.8 (10 orbital periods) 
except for the finer grid case, in which the simulation is 
run up to t = 12.56. We find that the main flow patterns 
reach a steady state by these times. 

3. Numerical results 

3.1. Astrostrophic wind 

Figures |l](a)-(d) show the equi-density surface of the com- 
panion star and the streamlines starting from the equi- 
density surface; Fig. |l](a) depicts the equi-density surface 
of log 10 p = —1, (b) —1.5, (c) —2, and (d) —2.5. In order 
to draw the streamlines, we use the same integration time 
step. Therefore, the higher the density of gas, the shorter 
the streamlines as the velocity is lower. We observe that 
there are ascending motions of gas in the higher density 
region (deeper inside the star) and a remarkable circu- 
lating flow in the lower density region (surface region). 
The streamlines of lighter color have higher velocity than 
those of darker ones. As will be shown later, there is a high 
pressure around the north pole, and so we may call this 
circulatory flow the north pole circulatory flow or simply 
an H-eddy. 

In order to see the surface flow more closely, Fig. |^ 
shows the equi-density surface (log 10 p = —2) of the com- 
panion star and the streamlines starting from the equi- 
density surface of log 10 p = —2.5 seen from various view 
angles. We confirm the presence of the H-eddy again. Gas 
in the H-eddy seems to drift gradually downward and 
its velocity increases with decreasing latitude. This phe- 
nomenon will be discussed later. 

Figure ^ shows the iso-velocity lines plotted on the 
equi-density surface of log 10 p = —2.5 seen from the north 
pole. The numbers in the figure show the magnitude of 
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Fig. 1. Equi-density surface of the companion star and streamlines starting from the equi-density surface. Integration 
time steps to calculate the streamlines are the same for all streamlines, and so the shorter the lines the slower the 
speed of gas. Note that the equi-density surface is not completely opaque. The equi-density surface are defined by (a) 
log 10 p = -l; (b) -1.5; (c) -2; (d) -2.5. 
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Fig. 2. Streamlines starting from the equi-density surface 
of log 10 p — — 2.5 and equi-density surfaces of log 10 p = — 2 
seen from various view angles. 




Fig. 3. Iso-velocity lines plotted on the equi-density sur- 
face of log 10 p = —2.5 seen from the north pole. The gray 
color shows the combined effect of the magnitude of the 
velocity and the lighting effect. 




Fig. 4. Isobaric lines plotted on the equi-potential surface 
seen from the north pole. 



the velocity (note that the velocity is normalized by Ail, 
where A is the separation and Q is the rotational frequency 
of the binary.) 

In the LI region, which is close to the LI point, there is 
a strong low pressure due to the mass outflow through the 
LI point. Because of the presence of the LI low pressure, 
there is a circulatory flow, around the LI point, rotat- 
ing counter-clockwise, although this eddy does not form a 
complete circle for geometrical reasons. We call this flow 
the Ll-eddy, although, strictly speaking, this is not a com- 
plete eddy. 

The other remarkable feature is the presence of two 
eddies rotating counter-clockwise at the opposite side of 



the LI point. These two eddies merge to form one larger 
eddy at different density levels. Since this region is close 
to the L2 point, we may call this eddy the L2 eddy. Since 
these H, LI and L2 eddies are explai ned i n terms of the 
astrostrophic wind (Lubow and Shu 1975 ), the pressure 
distribution on the surface of the companion is important. 

3.2. High / low pressures on the companion surface 

If there is no flow on the companion, the isobaric surface 
should coincide with the equi-potential surface as a re- 
sult of the hydrostatic balance. In our case, the isobaric 
surface does not coincide with the equi-potential surface 
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Fig. 5. Isobaric lines plotted on the equi-density surface 
of lo 

§io P — ^ 2.5 seen from the north pole. 



in an exact sense because of the presence of the surface 
flow. Fig. [| shows the isobaric lines plotted on the equi- 
potential surface of the companion. This equi-potential 
surface nearly corresponds to the equi-density surface with 
log 10 p = —2.5 except near the LI region, where high speed 
gas motion is present. Note that the isobaric lines on the 
companion surface have concentric patterns around the 
north pole. 

As was described earlier, the flow pattern on the sur- 
face of the rotating (in the inertial frame) companion can 
be explained in terms of the astrostrophic wind, which is 
a counterpart of the geostrophic wind in the atmosphere 
of the Earth. If viscosity can be neglected, the direction of 
the flow is parallel to the isobaric lines of low/high pres- 
sure. In the present case of a counter-clockwise rotating 
binary system (as viewed from the north pole), the flow 
proceeds by keeping the high (low) pressure to the right 
(left) in the northern hemisphere. Therefore, the H-eddy 
circulates clockwise, while the LI eddy and L2 eddy cir- 
culate in the counter-clockwise direction. 

If the pressure is only a function of the density of gas, 
the gas is called barotropic. Otherwise, the gas is called 
baroclinic. In the case of 7 = 1.01, the gas is nearly iso- 
thermal, hence almost barotropic. In the other cases, the 
gas is baroclinic. In a barotropic gas, an isobaric surface 
must coincide with an equi-density surface. Fig. || shows 
the baroclinicity of the gas with 7 — 5/3, whose property 
is responsible for the generation of vorticity. This will be 
discussed later. 



3.3. Ekman boundary layer 

In the Earth's atmosphere, an Ekman boundary layer ex- 
ists close to the surface of the Earth, where (turbulent) 
viscosity cannot be neglected and where the direction of 
the flow deviates from that of a pure geostrophic wind. In 



low pressure regions, gas drifts inwards and eventually as- 
cends near the centre of the low pressure regions (Pedlosky 
19871 ) . 

In our simulations, there is no (physical) viscosity as- 
sumed. Nevertheless, numerical viscosity arising from a 
truncation error acts as a kind of a turbulent viscosity, 
leading to gas drifting. This is the reason why the circulat- 
ing flow about the north pole gradually migrates towards 
the equatorial region. In order to test the above consider- 
ation, we perform a calculation with the mesh size halved. 
In this case the magnitude of the (numerical) viscosity is 
expected to be 1 /S of the coarser grid case. We find a simi- 
lar flow pattern as before, although the circulating motion 
is tighter than before (see Figs. ^). The main characteris- 
tics of the flow pattern are not changed, though. 

In the surface region of low-mass main-sequence stars, 
there is a convection zone, which provides the necessary 
turbulent viscosity. We may therefore expect that the 
present flow pattern is qualitatively correct. 

It is however not trivial from our simulations to de- 
rive the exact value of our (numerical) viscosity. We can 
only make a very tentative guess that the typical Reynolds 
number in our simulations is probably of the order of 10 to 
100. In a future work, we plan to solve this problem using 
the Navier-Stokes equation instead of the Euler ones, and 
using the molecular viscosity as a rough estimate of the 
turbulent viscosity which should be at play here. 

In our case the inner boundary does not play the role 
of the surface of the Earth. The Ekman layer is present at 
the top of the atmosphere, a situation similar to a water 
vortex. Gas circulating around the high pressure region 
at the north pole migrates gradually towards the equator, 
switches to the counter-clockwise motion close to the LI 
low pressure, and finally escapes from the companion to- 
wards the primary star. This flow pattern is important to 
understand the structure of the LI stream to be discussed 
later. 



3.4. Effect of 7 

Until now we have shown the results for 7 = 5/3, corre- 
sponding to ionized hydrogen gas. In order to study the 
case of molecular hydrogen, we calculate the case 7 = 7/5. 
The results are essentially the same as for 7 = 5/3, so that 
we do not discuss them any further. 

However, the results for 7 = 1.01 are very different 
from the cases of 7 = 5/3 and 7/5, and are essentially 
the same as that of Fujiwara et al. (2001): there is no H 
eddy but the flow originating at the pole region migrates 
rather directly towards the equatorial region. More pre- 
cisely, such a flow is rather weak, because the pressure on 
an equi-potential surface is rather uniform except in the 
region close to the LI point. The LI flow mainly originates 
from the region close to the LI point. The case of 7 = 1.2 
is intermediate to the above two extreme cases. We found 
the L2 eddy in this case. 
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Fig. 6. (a) Calculation with mesh size 121x101x51. (b) 
Calculation with mesh size 241x201x101 at t = 12.56. 



We may explain the difference of the flow patterns be- 
tween the cases of 7 = 5/3,7/5 and that of 7 = 1.01 in 
terms of the baroclinicity of the gas. Vorticity is gener- 
ated by the term containing Vp x V p, whi ch is identically 
zero for a barotropic gas (Pedlosky |l987|) . Therefore, lo- 
cal vorticity of the gas in the case of 7 = 1.01 is almost 
zero everywhere. Moreover, the Taylor-Proudman theo- 
rem, which states that the horizontal velocity components 
u and v are independent of z, holds not only for an incom- 
pressible flow but also for a barotropic gas (Busse 1970 ; 
Pedlosky |i~987l ). Since u and v are almost zero on the in- 
ner boundary, they must also be small in the upper at- 
mosphere. This problem, however, needs further investi- 
gation. We note that the cases of 7 = 5/3 and 7/5 are 
more appropriate to compare with (future) observations. 




Fig. 7. Iso-Mach number lines on an equi-density surface 
with log 10 p = —2.5 seen from the north pole. The sonic 
line is close to the LI point. The gray color shows the 
combined effect of the magnitude of the velocity and the 
lighting effect. The heavy black dot shows the position of 
the inner LI point. 



3.5. LI region and LI stream 



Lubow and Shu (1975) paid lots of effort to solve the flow 
in the LI region by their matched expansion technique. In 
this subsection we summarize our numerical results. 

In order to make the influence of the outer boundary 
to the LI stream small, we slightly expand the computa- 
tional region of the x-component from 0.7 to 0.9 in this 
subsection. The primary star is still out of the computa- 
tional region. Thus the LI stream reaches the boundary 
at x = 0.9 and escapes from the calculation region. In 
the present calculation no accretion disc is formed, and 
so there is no collision between the LI stream and the 
accretion disc aroun d the primary star. 

Lubow and Shu ( 1975 ) predicted that the sonic point, 
where subsonic gas is accelerated to supersonic velocity, 
must occur close to the LI point. Figure fj] shows the iso- 
Mach number lines on an equi-density surface and the 
numbers on the figure show the Mach number. It can be 
seen that the sonic line on the surface is certainly close to 
the LI point. 

Figs. ||(a)-(c) show the velocity vectors, the stream- 
lines, and the equi-density lines of the gas in the orbital 
plane, respectively. The flow pattern in the orbital plane 
near the LI point shows some difference with the one de- 
picted in Fig. 3 in Lubow and Shu ( 1975| ). They predicted 
that the streamlines of the lower portion of the LI point 
turn back to the companion, while in our case all gas flow 
towards the primary compact component. There are two 
differences between thier assumptions and ours: 1) they 
assumed iso-thermal gas; 2) they start their calculation 
at the LI point, where the gas is assumed to be at rest. 
The assumption 2) may be the main reason. We solve the 
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Fig. 8. (a) Velocity vector in the orbital plane. Note that 
arrows are drawn only every four cells in each x and y 
direction, (b) Streamlines in the orbital plane with veloc- 
ity magnitude contours, (c) Density contours with equi- 
density lines. The heavy black dot shows the position of 
the inner LI point. 



Fig. 9. Density contours of the LI stream at two cross sec- 
tions: (a) x = 0.5; (b) x = 0.65. The horizontal axis is the 
y-axis ranging from -0.5 to 0.5. The vertical axis depicts 
the z-axis ranging from to 0.5. Note the asymmetry of 
the density contour in the central part of the LI stream. 



whole region including both the LI region and the surface 
layer of the companion. The gas at the LI point is already 
supersonic as was described above. 

From Fig. ||(b) we can measure the deflection angle 9 S 
of the LI stream, which is the angle between the central 
flow line in the LI stream and the x— axis. It is about -10 
degrees. This angle was calculated to b e abou t -20 degrees 
for the case of q = 1 by Lubow and Shu ( 1975 ). In Fujiwara 
et al.'s result ( |200l| ) this is about -15 degrees. This differ- 
ence is also due to the difference in boundary conditions. 
Lubow and Shu assumed that the fluid element was at 
rest at the LI point. The effect of the Coriolis force is 
measured by the Rossby number defined by U/AQ, where 
U is a typical velocity of gas. The lower the initial velocity 
of the gas, the lower the Rossby number, the larger the 
Coriolis' effect and, hence, the deflection angle. 

In Fig. H(c) we found that the atmosphere of the sec- 
ondary is thicker in the leading hemisphere (negative y) 
than the trailing hemisphere. This can be easily explained 
by the direction of the Coriolis force, which acts to the 
right side of the flow direction in our system. 

Because of this asymmetry of the atmosphere, the 
structure of the LI stream is also asymmetric about the 
central flow line. In Figs. ^ we show the cross section of 
the LI stream, which is bow shaped convex to the positive 
y direction. We obtained a similar pattern in Fujiwara et 
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Fig. 10. Density profile along the y-axis on the orbital 
plane z — at various x locations, ranging from 0.5 to 
0.8. 



al. ( 2001 ), and this was explained by the interaction of 
the LI stream and the circulating circum-disc flow on the 
orbital plane about the primary. In the present work we 
realize that this explanation is not correct. The bow shape 
of the LI stream should be explained by the structure of 
the surface flow on the companion. 

Lubow and Shu (1975) predicted a symmetric 
Gaussian density distribution at the centre of the stream 
as 



p(r,d)^-exp[-Br 2 (6-6 s f 
r 



(1) 



where r is the distance from the LI point in the scaled co- 
ordinate, 9 S is the angle of the stream centre from the line 
joining the centre of the two stars, and C, B are the coeffi- 
cients depending on the mass ratio. Fig. flo| shows that our 
density profile nearly has a Gaussian density distribution. 
The density at the centre of the stream obtained by them 



decreases as 



while our result decreases as 



-0.8 



Fig. 



[lT| shows the equi-density surface of \ogp = —2.5 of the 
companion star and the LI stream. 

4. Discussion 

In our simulations, we obtained a system of eddies: the 
H eddy rotating clockwise, and counterclockwise rotating 
LI and L2 eddies. We now discuss what we believe is the 
reason that this eddy system is formed (see Fig. |l2| ). 

If we consider a time sequence, the LI low pressure 
is formed first, because of the mass loss through the LI 
point. Since the gas near the equator does not feel much 
resistance in flowing towards the LI point, the equato- 
rial region becomes a low pressure zone, and a high pres- 
sure zone builds up at the north pole: the H eddy is thus 
formed. 

The above explanation does not explain the appear- 
ance of the third eddy: the L2 eddy. This can be explained 
as follows. Consider a force exerted by an eddy on another 




Fig. 11. Equi-density surface of log 10 p = —2.5 of the com- 
panion star and the LI stream seen from the north pole. 
The heavy black dot shows the position of the inner LI 
point. 




Fig. 12. Schematic diagram of the system of eddies: LI, 
H, and L2 show the low pressure near the LI point, the 
high pressure around the north pole, and the low pressure 
at the opposite side of LI, respectively. An eddy rotates 
clockwise for H, and counterclockwise for LI and L2. 



eddy. The LI eddy tries to convect the H eddy to nega- 
tive y-direction, while the H eddy tries to convect the LI 
eddy to negative y-direction as well. However, the LI eddy 
must be fixed to the LI region, so that the H eddy must 
drift to negative y-direction. The system is therefore not 
steady. In order to stabilize the eddy system, we need the 
third eddy, which tries to convect the H eddy to positive 
y-direction. The L2 eddy thus appears and evolves until 
the three eddy system is stabilized. The exact location and 
the strength of these eddies can be obtained only through 
numerical simulations. 
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It is interesting to note that Davey and Smith (1992) 
have found that the trailing hemisphere has stronger line 
absorption than the leading hemisphere in a dwarf nova 
(see also Dhillon and Watson 2000). This asymmetry can 
be naturally explained in terms of the H-eddy current, 
which convects heat due to irradiation to the leading hemi- 
sphere. 

If this kind of astrostrophic wind conveys heat effec- 
tively from an irradiated hemisphere to the other, the unir- 
radiated side of the companion becomes hotter than pre- 
viously thought. Such an example may be observed in the 
2000 outburst of the recurrent nova CI Aql because, in 
nova outbursts, a companion star should be strongly irra- 
diated by a hot white dwarf. The inclination angle of the 
CI Aql orbit is about 80 degrees and the eclipse is almost 
total in the sense that the companion completely occults 
the bright accretion disc around the white dwarf (e.g., 
Hachisu & Kato 2001). In eclipse minima, therefore, we 
see only the unirradiated side of the companion. Hachisu, 
Kato and Schaefer (2002) showed that the eclipse min- 
ima of CI Aql in outburst is about 0.6 mag brighter than 
in quiescence and could reasonably reproduce the orbital 
light curve (folded by the orbital period) when the tem- 
perature of the unirradiated hemisphere is about 1000 K 
hotter than that in quiescence. 

In the present study, we do not take into account such 
complex effects as (real) viscosity, magnetic fields, and ir- 
radiation from the accretion disc. Activities of magnetic 
fields on the companion surface would lead to star-spots, 
and the irradiation from the accretion disc would lead 
to an increase of the gas temperature of the companion. 
These effects might more or less change the surface flow 
patterns. The present study is a first step to examine the 
surface flow structure. Thus to study these phenomena 
is beyond the scope of the present paper, and is left for 
future work. 



5. Summary 

We performed numerical simulations of the surface flow 
on the companion star in a semi-detached binary system. 
The case with 7 = 5/3 was mainly studied, although the 
cases of 7 = 7/5, 1.2 and 1.01 were also investigated. The 
results are summarized as follows: 

1. We obtain an eddy configuration composed of an H- 
eddy, an Ll-eddy, and an L2-eddy associated with 
high/low pressures. If there is no viscosity, the flow 
lines must be closed. The above eddies are explained 
in terms of the astrostrophic wind. 

2. Some kind of viscosity (e.g., turbulent viscosity origi- 
nating from a surface convection zone) would lead to 
a downward drifting motion. In our simulation the nu- 
merical viscosity acts as turbulent viscosity. 

3. The sonic line on an equi-density surface occurs close 
to the LI point. 



4. The deflection angle of the LI stream is about -10 to 
-15 degrees. This small angle is due to the larger fluid 
velocity at the LI point. 

5. The density distribution in the LI stream is nearly 
symmetric about the central flow line. The cross sec- 
tion of the LI stream is bow shaped convex to the 
positive y direction. This is an intrinsic nature of the 
stream and is not caused by the interaction with the 
accretion disc gas. 
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